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Abstract 

Geo-neutrino detection can determine the amount of long-lived radioactive elements within our 
planet, thus providing a direct test of the Bulk Silicate Earth (BSE) model and fixing the radio- 
genic contribution to the terrestrial heat. We present a prediction for the geo-neutrino signal at 
KamLAND as a function of the Uranium mass in the Earth. The prediction is based on global 
mass balance, supplemented by a detailed geochemical and geophysical study of the region near 
the detector. The prediction is weakly dependent on mantle modeling. If BSE is correct, Uranium 
geo-neutrinos will produce between 25 and 35 events per year and 10 32 protons at Kamioka. 
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I. INTRODUCTION 



The deepest hole that has ever been dug is about 12 Km deep, a mere dent in planetary- 
terms. Geochemists analyze samples from the Earth's crust and from the top of the mantle, 
whereas seismology can reconstruct the density profile throughout the whole Earth but not 
its composition. In this respect, our planet is mainly unexplored. Geo- neutrinos — the 
antineutrinos from the progenies of U, Th and 40 K decays in the Earth — bring to the 
surface information from the whole planet, concerning its content of radioactive elements. 
Their detection can shed light on the sources of the terrestrial heat flow, on the present 
composition and on the origins of the Earth. 

Geo-neutrino properties, summarized in Table EJ deserve a few comments: 

(i) geo-neutrinos originating from different elements can be distinguished due to their 
different energy spectra, e.g., geo-neutrinos with energy E > 2.25 MeV are produced only 
from the Uranium chain; 

(ii) geo-neutrinos from U and Th (not those from 40 K) are above threshold for the classical 
anti-neutrino detection reaction, the inverse beta on free protons: 

u + p -»• e + +n- 1.8 MeV ; (1) 

(iii) anti-neutrinos from the Earth are not obscured by solar neutrinos, which cannot 
yield reaction (JTJ). 

TABLE I: The main properties of geo-neutrinos: en is the heat production rate per unit mass 
and natural isotopic composition; ep is the antineutrino production rate (number of antineutrinos 
per unit time) per unit mass (Contribution of 235 TJ is neglected due to the small, 0.7%, natural 
abundance). 

Decay Q r 1/2 E max 

[MeV] [10 9 yr] [MeV] [W/Kg] [kg^s" 1 ] 
238tj _> 206p b + 8 4 He + 6e + 6f? 51 j 4A7 3 2 6 0.95 x 10~ 4 7.41 x 10 7 
232 Tll _ > 208 Pb + 6 4 He + 4e + 4z7 42 7 140 2.25 0.27 x 1(T 4 1.63 x 10 7 
40 K _^ 40 Ca + e + p 1.32 1.28 1.31 0.36 x 1(T 8 2.69 x 10 4 
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Geo-neutrinos were introduced by Eder (lj in the sixties and Marx soon realized their 
relevance. In the eighties Krauss et al. discussed their potential as probes of the Earth's 
interior in an extensive publication Q]. In the nineties the first paper on a geophysical 
journal was published by Kobayashi et al. [4]. In 1998, Raghavan et al. jjj and Rotschild et 
al. (f| pointed out the potential of KamLAND and Borexino for geo-neutrino detection. 
In the last three years more papers appeared than in the previous decades: in a series of 



papers Fiorcntini ot al. [tUs!,!^ 10 




11 1 discussed the role of geo-neutrinos for determining the 



radiogenic contribution to the terrestrial heat flow and for discriminating among different 
models of Earth's composition and origin. A reference model for geo-neutrino production, 
based on a compositional map of the Earth's crust and on geochemical modeling of the 
mantle, was presented in At the end of 2002, the analysis of the first data release of 
KamLAND |l2j (equivalent to an exposure of 0.11 x 10 32 proton yr and 100% efficiency) 
reported 4 events from Uranium and 5 from Thorium out of a total of 32 counts in the 
geo-neutrino energy region (E v i s < 2.6 MeV), after subtracting 20 reactor events and 3 
background counts. Statistical fluctuations imply an error of, at least, 5.7 counts. Indeed 



this first indication of geo-neutrinos stimulated several investigations [13 . 



191, 

In a few years KamLAND should provide definite evidence of geo-neutrino signal, after 
accumulating a much larger statistics and reducing background. In the meanwhile other 
projects for geo-neutrino detection are being developed. Borexino at Gran Sasso, which 
is expected to take data in a few years, will benefit from the absence of nearby reactors. 



Domogatski et al. 



2 it ] are proposing a 1 Kton scintillator detector in Baksan, again very 



far from nuclear reactors. A group at the Sudbury Neutrino Observatory in Canada is 
studying the possibility of using liquid scintillator after the physics program with heavy 
water is completed. The LENA proposal envisages a 30 Kton liquid scintillator detector 
at the Center for Underground Physics in the Pyhasalmi mine (Finland). Due to the huge 
mass, it should collect several hundreds of events per year. In conclusion, one can expect 
that within 10 years the geo-neutrino signal from Uranium and Thorium will be measured 
at a few points on the globe. 

In this paper we shall concentrate on geo-neutrinos from Uranium, which are closer 
to experimental detection, and on the predictions for Kamioka, the site hosting the only 
detector which is presently operational. Our goal is to understand which information on the 
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total amount of Uranium in the Earth can be extracted from geo-neutrino measurements. 

As briefly discussed in the next section, the Uranium mass in the Earth is estimated on 
the grounds of cosmo-chemical arguments, based on the compositional similarity between 
Earth and carbonaceous chondrites. Measurements of samples from the Earth's crust imply 
that the crust contains about one half of this global estimate, whereas the mantle — which 
should contain the rest — is practically unexplored in this respect. A direct determination 
of the Uranium mass in the globe is clearly an important test of the origins of the Earth. 
Furthermore, such a determination will also fix the radiogenic contribution to the terrestrial 
heat flow, which is a presently debated issue, see, e.g., Ref. j^]. Early estimates of the 
geo-neutrino signal and their connection with the global Uranium content are also reviewed 
at the end of Section HTl 

The geo-neutrino signal depends on the total Uranium mass m in the Earth and on the 
geochemical and geophysical properties of the region around the detector For Kam- 
LAND, we estimated [2| that about one half of the signal is originated within 200 km from 
the detector. This region, although containing a globally negligible amount of uranium, 
produces a large contribution to the signal as a consequence of its proximity to the detec- 
tor. This contribution has to be determined on the grounds of a detailed geochemical and 
geophysical study of the region, if one wants to extract from the total signal the remaining 
part which carries the relevant information on m. The study of the region around Kamioka, 
including the possible effects of the subducting plates below the Japan Arc and a discussion 
of the contribution from of the Japan Sea, is presented in Section IIIII 

The contribution from the rest of the world, discussed in Section IIV| depends on the 
total mass of Uranium as well as on its distribution inside the Earth, since the closer is 
the source to the detector the larger is its contribution to the signal. For each value of 
m, we shall construct the distributions of Uranium abundances which provide the maximal 
and minimal signals, under the condition that they are consistent with geochemical and 
geophysical information on the globe. For the Earth's crust we shall use a 2° x 2° map [^J 
which distinguishes several crustal layers and to each layer we shall assign minimal/maximal 
values for the Uranium mass abundances. According to geochemists, the rest of Uranium 
should be found in the mantle. Observational data for this reservoir are very poor, however 
it is generally believed that Uranium abundance increases with depth. The assumption 
that abundance is spherically symmetrical and non decreasing with depth will be enough to 
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provide rather tight bounds on the mantle contribution to the geo-neutrino signal. This will 
be further combined with the results for the crust. 

We shall put together the pieces of the above analysis in Section |VJ where we present 
our main result in Fig. ^ a narrow band describes the predicted signal as a function of 
Earth's total Uranium mass. We remark that the extremes of the band correspond to the 
whole range of uncertainty, which is estimated according to the following criteria: (i) for 
statistical errors we consider a ±3<r interval; (ii) for systematic uncertainties of geochemical 
and geophysical parameters we determine an interval such as to cover all modern estimates 
which we found in the literature; (iii) independent errors are combined in quadrature. 

On the grounds of Fig. ^ we discuss how the geo-neutrino signal can provide a direct test 
of a fundamental paradigm on the origins and composition of our planet. 

H R (U) [TW] 

5 10 15 




0.5 1 1.5 

m [10 17 kg] 



FIG. 1: The predicted signal from Uranium geo-neutrinos at KamLAND. 
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II. U, TH AND K IN THE EARTH: HOW MUCH AND WHERE? 



Earth's global composition is generally estimated from that of chondritic meteorites by 
using geochemical arguments which account for losses and fractionation during planet for- 
mation. Along these lines the Bulk Silicate Earth (BSE) model is built, which describes 
the "primitive mantle", i.e., the outer portion of the Earth after core separation and before 
the differentiation between crust and mantle. The model is believed to describe the present 
crust plus mantle system. Since lithophile elements should be absent in the core 1 , the BSE 
provides the total amounts of U, Th and K in the Earth, estimates from different authors 
being concordant within 10-15% [2^. From the estimated masses, the present radiogenic 
heat production rate H R and anti-neutrino luminosity L v can be immediately calculated, 
see Table HTl and, e.g., Ref. jlpl |. 



TABLE II: U, Th and K according to BSE, from Ref. 
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The BSE is a fundamental geochemical paradigm. It is consistent with most observations, 
which however regard mostly the crust and an undetermined portion of the mantle. The 
measurement of quantities — such as the geo-neutrino signals — which are directly related 
to the global amounts of radioactive elements in the Earth will provide a direct test of this 
model for composition and origin of our planet. 

For sure, heat released from radiogenic elements is a major source of the terrestrial heat 
flow, however its role is not understood at a quantitative level. The masses estimated 
within the BSE account for the present radiogenic production of 19 TW, which is about one 

1 One needs to be careful, since the definition of an element's behavior, i.e., lithophile or not, depends on 
the surrounding system; there exist models of the Earth's core suggesting it as repository for radioactive 
elements. 
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half of the estimated heat flow from Earth 23J, |24| . Anderson refers to this difference as the 



missing heat source mystery and summarizes the situation with the following words: "Global 
heat flow estimates range from 30 to 44 TW . . . Estimates of the radiogenic contribution 
. . .based on cosmochemical considerations, vary from 19 to 31 TW. Thus, there is either a 
good balance between current input and output ... or there is a serious missing heat source 
problem, up to a deficit of 25 TW ..." If one can determine the amounts of radioactive 
elements by means of geo-neutrinos, an important ingredient of Earth's energetics will be 
fixed. 

Concerning the distribution of radiogenic elements, estimates for Uranium in the conti- 
nental crust based on observational data are in the range: 

m c = (0.3 - 0.4) x 10 17 kg . (2) 

The extreme values have been obtained in Ref. |lj] by taking the lowest (highest) concen- 
tration reported in the literature for each layer of the Earth's crust, see Table II of Ref. jsj], 
and integrating over a 2° x 2° crust map. The main uncertainty is from the Uranium mass 
abundance a^c in the lower crust, with estimates in the range (0.2 — 1.1) ppm. Estimates 
for the abundance in the upper crust, due, are more concordant, ranging from 2.2 ppm to 
2.8 ppm. The crust — really a tiny envelope — should thus contain about one half of the 
BSE prediction of Uranium in the Earth. 

About the mantle, observational data are scarce and restricted to the uppermost part, 
so the best estimate for its Uranium content uim is obtained by subtracting the crust con- 
tribution from the BSE estimate: 

m M = m B sE - m c . (3) 

A commonly held view is that there is a vertical gradient in the abundances of incompat- 
ible elements in the mantle, with the top being most depleted. A minimum gradient model 
has a fully mixed and globally homogeneous mantle; the other extreme is a model where all 
the Uranium is at the bottom of the mantle. 

Geochemical arguments are against the presence of radioactive elements in the completely 
unexplored core, as discussed by McDonough in a recent review of compositional models of 
the Earth 



S 



Similar considerations hold for Thorium and Potassium, the relative mass abundance 
with respect to Uranium being globally estimated as: 



a(Th) : a(U) : a(K) « 4 : 1 : 10, 000 . (4) 

We remark that the well-fixed ratios in Eq. (|IJ) imply that detection of geo-neutrinos from 
Uranium will also bring important information on the amount of Thorium and Potassium 
in the whole Earth. 

Several predictions for the geo-neutrino signal have been presented in the past, corre- 
sponding to different hypotheses about the amount of Uranium in the Earth and to different 
models of its distribution. A summary is presented in Fig. |21 Early models (full 
circles) assumed a uniform Uranium distribution in the Earth and different values of the 
Uranium mass. In fact these predictions are almost proportional to m. The huge signals 
predicted by Eder and by Marx were obtained by assuming that the Uranium density in 
the whole Earth is about the same as that observed in the continental crust; Marx (Eder 2 ) 
assumed thus an Uranium mass 30 (60) times larger than that estimated within the BSE. 

Krauss et al. (^J distributed about 10 17 kg of Uranium uniformly over a 30 km crust. 
The other estimates (crosses) are all obtained by using the BSE value for the mass as an 
input and different models for distributing the Uranium content between crust and mantle. 
In this class, Rotschild et al. Jf| obtained the minimal prediction by assuming for the crust 
a very small Uranium abundance, definitely lower than the values reported in more recent 
and detailed estimates. 

In this paper we shall use a rather general approach, by keeping the total amount of 
Uranium as a free variable, within the loose constraints provided from one side by the 
amount observed in the crust and from the other side by the amount tolerated by Earth's 
energetics. We shall distribute the total amount between crust and mantle so as to maximize 
or minimize the signal, within the boundary provided by geochemical and/or geophysical 
observations. 



2 The factor 10.7 appearing in the last of the equations (13) of Ref. [l| should actually be 1.07 and, 
correspondingly, the reaction rate on page 661 has been divided by a factor of 10. 
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FIG. 2: Previous estimates of the geo-neutrino signal S, renormalized to the average survival 
probability (P ee ) = 0.59, and the corresponding estimated Uranium mass m. The signal is in 
Terrestrial Neutrino Units (1 TNU = 1 event /year /10 32 proton). 

III. THE REGION NEAR THE DETECTOR 



As mentioned in the introduction, the entire Earth's crust will be subdivided into 2° x 2° 
tiles. Within each tile, one distinguishes several vertical layers and assigns to each layer a 
world averaged Uranium mass abundance, see Ref. [jj. With the aim of reducing the error 
on the regional contribution to the level of the uncertainty on the rest of the world, one 
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needs a more detailed geochemical and geophysical study of the crust in the region within a 
few hundreds kilometers from the detector, where some half of the signal is generated. 

We shall present here our results for the region near the KamLAND detector, located at 
36° 25' 26" N and 137° 19' 11" E. We analyze the six tiles (see Fig.EJ) around KamLAND by 
using geochemical information on a 1/4° x 1/4° grid and a detailed map of the crust depth. 
The possible (minimal and maximal) effects of the subducting slab beneath Japan are also 
considered and the uncertainty arising from the debated (continental or oceanic) nature of 
the crust below the Japan Sea is taken into account. 



A. The six tiles near KamLAND 



The seismic velocity structure of the crust beneath the Japan Islands has been determined 
in Ref. from the study of some 13,000 arrival times of 562 local shallow earthquakes. By 
applying an inversion method, the depth distribution of the Conrad and Moho discontinuities 
beneath the whole of the Japan Islands are derived, with an estimated standard error of 
±1 km over most of Japan territory. Our Figs. El and 0] are derived from Fig. 6 of Ref. (27 1. 
This allows distinguishing two layers in the crust: an upper crust extending down to the 
Conrad and a lower part down to the Moho discontinuity. 

The upper-crust chemical composition of Japan Islands has been studied in Ref. (^J, 
based on 166 representative specimens, which can be associated with 37 geological groups 
based on ages, lithologies and provinces. By combining the base geological map of Fig. 2 
of Ref. j3] — which distinguishes 10 geological classes — with the abundances reported in 
Table 1 of the same paper, one can build a map of Uranium abundance in the upper crust, 
under the important assumption that the composition of the whole upper crust is the same 
as that inferred in Ref. |28| from the study of the exposed portion, see Fig. |SJ 

We are not aware of a specific study of the lower part of the Japan crust, however, it is 
well known that there are similarities between the composition of the Japanese crust and 
that of the Sino-Korean block. In an extensive compositional study of East China crust J^], 
the Uranium abundance in the lower part is estimated between 0.63 and 1.08 ppm. On these 
grounds we shall take for the abundance in the lower crust of Japan: 

a LC = 0.85 ±0.23 . (5) 
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FIG. 3: Depth of the Conrad discontinuity in Japan, from Ref. 
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We remark that the estimated range of abundances for East China is substantially nar- 
rower than the range for the whole world, which is (0.2-1.1) ppm (see Table II in Ref. j^j]), 
this last interval presumably reflecting regional differences in the lower crust composition. 

Concerning the vertical distribution of abundances in the crust, it is presently impossible 
to have information on the chemical composition on a scale smaller than the Conrad depth, 
generally lying at about 20 km. We also note that the sampling density for the study 
of the upper crust in the region near Kamioka is about one specimen per 400 km 2 . On 
these grounds, we introduce a grid where one assigns a specific abundance to cells with size 
1/4° x 1/4°, i.e., with a linear scale of about 20 km. Within each cell, the depth of the upper 
and lower crust are taken from Ref. 22J, Uranium abundance for the upper crust is derived 
from Ref. 



28| and for the lower crust from Ref. 29]. In this way each of the six tiles near 
Kamioka is subdivided into sixty-four cells. 

Just for the sake of computing flux and signal, each cell is further subdivided into many 
subcells with abundance derived from those of the parent cell. The angle-integrated produced 
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FIG. 4: Depth of the Moho discontinuity in Japan, from Ref. 



flux at distance R from a subcell of volume AV, calculated as that from a sphere with radius 
r = (3/47TAV) 1 / 3 , is 

\R-r\~ 



A$ = — 
AR 



2Rr + (R 2 -r 2 )ln- 



R + r • < 6) 
where A is the specific activity (number of neutrinos produced per unit time and volume). 
Each subcell provides a contribution to the signal rate AS" given by: 



AS = N p A$(aP ee (R)) 



(7) 



where N p is the number of free protons in the target, a is the cross section of reaction in 
Eq. |T|). P ee (R) is the survival probability which we shall calculate for tan 2 9 = 0.40 and 
Am 2 = 7.9 x 10 -5 eV 2 [3^. The average is over the energy spectrum of the neutrinos from 
Uranium decay chain. 

The resulting signal is obtained by adding the contributions of all subcells and will be 
expressed in Terrestrial Neutrino Units (TNU), where 1 TNU corresponds to 1 event per 
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FIG. 5: Uranium abundance in the upper crust of Japan 

year and per 10 32 protons. The contributions to the produced flux and to the signal from 
the six tiles are: 

$ 6 = 1.59 x 10 6 cm-V 1 ; (8a) 
S 6 = 12.74 TNU . (8b) 

With respect to our previous estimate from the whole globe [9j, giving <!> = 3.676 x 
10 6 cm -2 s _1 and S = 28.6 TNU for (P ee ) = 0.59, we find that the six tiles contribute 
43% of the flux and 45% of the signal: this justifies the close scrutiny of the region within 
the six tiles. Some 3/4 of the contribution arises form the upper crust. 

In more detail, the tile hosting Kamioka generates 29% and 30% of the total produced flux 
and signal, respectively. The host cell, i.e., the cell where Kamioka is located, contributes 
9% to the total produced flux. 

The Uranium mass contained in the six tiles is about ttlq = 3.3 x 10 13 kg, really negligible 
(less then 0.05%) with respect to that estimated for the whole Earth. 

We have considered several sources of the uncertainties affecting this estimate, see Ta- 
ble nn 
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(i) Measurement errors of the chemical analysis. The standard error of the Uranium 
abundance measurements in the individual samples is 3-4% [28]. We translated this error 
into a ±10% global uncertainty (correspondingly to about 3a) on the Uranium abundance 
in the upper crust for the six tiles 3 . 

(ii) Discretization of the upper crust. As mentioned above, we divided the crust into 
1/4° x 1/4° cells, assigning specific abundances to each cell. This discretization procedure 
introduces some uncertainty, which is especially important in the region very close to the 
detector. We have evaluated the effect of replacing the abundance of the host cell with those 
of adjacent cells. This produces signal variations in the interval (—0.64, +1.68) TNU. For 
simplicity, we introduce a symmetrical error such as to encompass the extreme values. 

(Hi) Chemical composition of the lower crust. The error is taken as the half- difference 
between signals obtained for the extreme values of the estimated Uranium abundances in 
the lower crust. 

(iv) Crustal depth. Since the depths of the Conrad and the Moho discontinuities are 
estimated with an 3a accuracy of about ±3 km, we have evaluated the effect of such (global) 
variations over the six tiles, the error being again estimated as the half difference between 
extreme values of flux/signal. 

TABLE III: Errors from the regional geophysical and geochemical uncertainties. 



Source 


AS (TNU) Remarks 


Composition of upper-crust samples 


0.96 


3a error 


Upper-crust discretizaion 


1.68 




Lower-crust composition 


0.82 


Full range 


Crustal depths 


0.72 


3a error 


Subducting slab 


2.10 


Full range 


Japan Sea 


0.31 


Full Range 


Total 


3.07 





3 This choice is very conservative given our lack of information on the correlation between the errors of 
the 166 samples: the errors should partially average out and result in a total error between 12% and 
(12/ v / 166)% « 1%. 
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B. The effect of the subducting slab beneath Japan 



The Japan arc, at the crossing among the Eurasian, Philippine and Pacific plates, is 
the theater of important subduction processes. The Philippine plate is moving towards the 
Eurasia plate at about 40 mm/yr and is subducting beneath the southern part of Japan. The 
Pacific Plate is moving in roughly the same direction at about 80 mm/yr and is subducting 
beneath the northern half of Japan. 

We shall model these processes as a single slab penetrating below Japan with velocity 
v = 60 mm/yr, the average of the two plates. This process has been occurring on a time 
scale T pa 10 s y, during which the front has advanced by D = vT pa 6000 km. 

We assume that the slab brings with it a sediments layer (with density p scd pa 1.6 ton/m 3 , 



depth h se d ~ 35C 



m and Uranium abundance a sed = 1.4 ppm, according to the data for the 



Japan trench pL\) on top of an oceanic crust layer, with density poc ~ 2.9 ton/m 3 , vertical 
extension hoc ~ 6.5 km, and Uranium abundance ape — 0.1 ppm. 
The amount of Uranium carried by the slab per unit surface is thus: 

cr s i ah = a scd x p sed x h sed + a oc x p oc x h oc = 2.7 kg/m 2 ■ (9) 

We observe that the corresponding quantity for the lower continental crust of Japan 
(density plc ~ 2.7 ton/m 3 , average depth Hlc ~ 19 km, and Uranium abundance a^c = 
0.85 ppm) is: 

olc = a LC x p LC x h LC = 43.6 kg/m 2 . (10) 

In order to estimate the effect of the subducting slab on geo-neutrino production, one 
can envisage two extreme cases: 

(a) one assumes that the slab keeps its trace elements while subducting. The effect of 
its presence can be estimated as if the lower crust is effectively enriched by the amount of 
Uranium contained in the crustal part of the slab passing below it, i.e., aic — > &lc + 0"siat» 
a negligible effect in comparison with the 25% uncertainty on a^c resulting from Aa^c ~ 
0.2 ppm. The signal is increased by 0.2 TNU. 

(b) At the other extreme, it is possible that, as the slab advances, all Uranium from 
the subducting crust is dissolved in fluids during dehydration reactions and accumulates in 
the lower part of the continental crust of Japan, thus strongly enriching it. The release 
process — to a first approximation — will be uniform along the subduction direction, for 
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some distance d ~ 250 km corresponding to the dimension of the Japan arc transverse to 
the trench. Since the slab has advanced by D « 6000 km, the abundance increase in the 
lower crust is now: 

Svlc = ^siab x ^ » 68.4 kg/m 2 . (11) 

This corresponds to a substantial increase of the effective Uranium abundance in the 
Japanese lower continental crust. The signal contributed from the lower crust of Japan 
increases by 4.4 TNU. 

As we have no argument for deciding which of the extreme cases (a) or (b) is closer 
to reality and in order to encompass both of them, we estimate the contribution from the 
subducting slab as: 

SWab — (2-3 ± 2.1) TNU . (12) 



C. The crust below the Japan Sea 

The morphology of the Japan Sea, characterized by three major basins (Japan, Yamato, 
and Ulleung Basins) and topographic highs such as the Yamato Ridge, is suggestive of 
intricate back-arc opening tectonics. Based on seismic reflection/refraction survey data, 
bottom sampling data, geomagnetic data and basement depth and topography, Tamaki et 
al. Q distinguish four crustal types: continental, rifted continental, extended continental, 
and oceanic crust. The Japan basin is generally considered as oceanic, whereas the nature 
of other basins is controversial and debated. Again, we resort to two extreme models: 

(a) following Ref. 2^ we consider all the basins as formed with oceanic crust, extending 
down to 7 km below 1 km of sediments. This provides a model for minimal geo-neutrino 
production, resulting in Sjs = 0.06 TNU. 

(b) Deeper crustal depths (up to 19 km for the Oki bank) and thicker sediments layers (up 
to 4 km for the Ulleung basin) are reported in the literature, see Table IIVI By taking these 
values and assigning the abundances typical of continental crust, we maximize geo-neutrino 
production, with the result Sjs = 0.68 TNU. 

In between (a) and (b), and in order to encompass the extreme values, we fix the contri- 
bution to the signal from the Japan Sea as: 

Sjs = (0.37 ± 0.31) TNU . (13) 
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TABLE IV: The vertical extensions (km) of crustal layers in the Yamano basin (YB), Oki bank 
(OK), and Ulleung basin (UB) used for model (b). 

YB OK UB 
sediments 1.2 0.3 4 
upper 2.8 8.7 2 
lower 8.5 10.5 8 

D. Summary of the regional contribution 

The regional contribution to the signal can be determined by adding the previous results. 
As the errors are independent, we combine them in quadrature, obtaining: 

S reg = (15.41 ± 3.07) TNU . (14) 

The principal uncertainty comes from the effect of the subducting slab. A more detailed 
study of the mechanisms of Uranium release should exclude the extreme cases which we 
have considered, thus reducing the error. 

Discretization of the crust is major uncertainty. A more detailed description of the 
exposed crust is certainly achievable, however, it will bring little help without a better 
understanding of the chemical composition variation with depth. 

IV. THE REST OF THE WORLD 

The contribution from the rest of the world will depend on the total mass of Uranium m 
as well as on its distribution inside the Earth, since the closer is the source to the detector 
the larger is its contribution to the signal. For each value of m, we shall construct the 
distributions which provide the maximal and minimal signals under the condition that they 
are consistent with geochemical and geophysical information on the globe. 

For the Earth's crust, we use the 2° x 2° map of Ref. Q distinguishing several crustal 
layers which are known to contain different amounts of radioactive elements. For each layer 
minimal and maximal estimates of Uranium abundances found in the literature are adopted, 
so as to obtain a range of acceptable fluxes, see Table 
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TABLE V: Minimal and maximal estimated Uranium abundances for the continental crust, in 
ppm. 





min 


max 


upper crust 


2.2 


2.8 


lower crust 


0.2 


1.1 



Depending on the adopted values, the Uranium mass in the crust mc is the range (0.3 — 
0.4) in units — here and in the following — of 10 17 kg. Clearly the larger is the mass the 
bigger is the signal, the extreme values being 4 : 

Sg° in) = 6.448 form = 0.3 and ^ max) = 8.652 form = 0.4. (15) 

Concerning Uranium in the mantle, we assume that spherical symmetry holds and that 
the Uranium mass abundance is a non- decreasing function of depth. It follows that, for a 
fixed Uranium mass in the mantle tum, the extreme predictions for the signal are obtained 
by: (i) placing Uranium in a thin layer at the bottom and (ii) distributing it with uniform 
abundance over the mantle. These two cases give, respectively: 

S { ™ n) = 12.15 x m M TNU and S { ™ x) = 17.37 x m M TNU . (16) 

By using again the proximity argument, we can combine the contributions from crust 
and mantle so as to obtain extreme predictions: for a fixed total m = mc + rriM, the highest 
signal is obtained by assigning to the crust as much material as consistent with observational 
data (mc = 0.4) and putting the rest, m — mc, in the mantle with a uniform distribution. 
Similarly, the minimal flux/signal is obtained for the minimal mass in the crust (mc = 0.3) 
and the rest in a thin layer at the bottom of the mantle. In conclusion, the contribution 
from the rest of the world is in the range: 

= [6.448 + 12.15(m - 0.3)] TNU and S { ™ x) = [8.652 + 17.37(m - 0.4)] TNU . 

(17) 

4 As we are now considering distances from the detector which are considerably larger than the neutrino 
oscillation length, the asymptotic expression for the survival probability, (Pee) = (1 — l/2sin 2 (2#)) = 
(1 + tan 4 6)/(l+ tan 2 0) 2 holds, so that the produced flux and signal are directly proportional, S/TNU = 
13.2 x (P ee )$/(f0 6 cm- 2 s- 1 ) Q. 
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Note that the two straight lines cross near m = 0.21, i.e., in the nonphysical region, since 
Uranium mass is at least 0.3. 

We remind that that both the Uranium poor (mc = 0.3) and the Uranium rich (mc = 0.4) 
crust models are observationally acceptable. We also recall that a compositionally uniform 
mantle is advocated by geophysicists, whereas geochemists prefer a two-layered mantle with 
a lower part close to the primitive composition and an upper part strongly impoverished in 
Uranium. In other words, it seems to us that the extreme predictions correspond to equally 
plausible models. On these grounds, we take as our prediction the mean of the extremes 
and assign an error so as to encompass both of them: 

S RW = [(2.25 + 14.76 x m) ± (-0.55 + 2.61 x m)] TNU . (18) 

We remark that by combining global mass balance with geometry, we have strongly 
constrained the contribution from the rest of the world: for a mass near the BSE values, 
m « 0.8, the signal is predicted within about 10%. 

V. THE GEO-NEUTRINO SIGNAL AS A FUNCTION OF URANIUM MASS IN 
THE EARTH 

By combining the regional contribution, Eq. (fPfj) . with that from the rest of the world, 
Eq. (|TB*|) . we get the Uranium geo-neutrino signal as a function of Uranium mass in the 
Earth: 

S = S ± A , (19a) 

where: 

S = 17.66 + 14.76 x m (19b) 
A 2 = (3.07) 2 + (2.61 x m-0.55) 2 . (19c) 

This error is obtained by combining in quadrature all geochemical and geophysical un- 
certainties discussed in the preceding paragraphs. All of them have been estimated so as 
to cover ±3cr intervals of experimental measurements and total ranges of theoretical predic- 
tions. 
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However, this error does not account for present uncertainties on neutrino oscillation 
parameters and on the cross section of reaction in Eq. (0). For the sake of discussing the 
potential of geo-neutrinos, we shall ignore for the moment these error sources. 

The expected signal from Uranium geo-neutrinos at KamLAND is presented as a function 
of the total Uranium mass m in Fig. ^ The upper horizontal scale indicates the correspond- 
ing radiogenic heat production rate from Uranium (Hr = 9.5 x m). 

The predicted signal as a function of m is between the two lines denoted as Shi g h and 
Slow, which correspond, respectively, to S ± A. 

Since the minimal amount of Uranium in the Earth is 0.3 x 10 17 kg (corresponding to 
the minimal estimate for the crust and the assumption of negligible amount in the mantle), 
we expect a signal of at least 19 TNU. On the other hand, the maximal amount of Uranium 
tolerated by Earth's energetics 5 , 1.8 x 10 17 kg, implies a signal not exceeding 49 TNU. 

For the central value of the BSE model, m = 0.8 x 10 17 kg, we predict S = 29.5±3.4 TNU, 
i.e., with an accuracy of 12% at "3cr". We remark that estimates by different authors for 
the Uranium mass within the BSE are all between (0.7 — 0.9) x 10 17 kg. This implies 
that the Uranium signal has to be in the interval (24.7 — 34.5) TNU. The measurement of 
geo-neutrinos can thus provide a direct test of an important geochemical paradigm. 

The effect of uncertainties about the oscillation parameters is presented in Table IVII In 
this respect the mixing angle is most important. Figure 4 (b) of Ref. [3^ shows a 3cr range 
0.26 < tan 2 9 < 0.67 (central value 0.40): the corresponding range for the average survival 
probability is 0.52 < P ee < 0.67 (central value 0.59), with a 3a relative error on the signal 
AS/S ~ 13%, which is comparable to the geological uncertainty in Eq. fJ19c|) . The mixing 
angle should be determined more precisely for fully exploiting the geo-neutrino signal. 

On the other hand, the predicted signal is practically unaffected by the uncertainty on 
the neutrinos squared mass difference Sm 2 : when this is varied within its ±3cr interval the 
signal changes by 0.1 TNU. This holds for any value of the total Uranium mass m, since the 
precise value of Sm 2 only matters in the region near the detector. In addition, we observe 
that the predictions computed for the best value (8m 2 = 7.9 x 10~ 5 eV 2 ) and for the limit 
Sm 2 = 00 differ by +0.3 TNU. Finally, the error on the inverse-beta cross section (quoted 

5 For an Uranium mass m = 1.8 x 10 17 kg and relative abundances as in Eq. the present radiogenic 
heat production rate from U, Th and K decays equals the maximal estimate for the present heat flow 
from Earth, iJ|? ax = 44 TW l^. 
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TABLE VI: Effect of the oscillation parameters on the signal. The relative/absolute variation 
is computed with respect to the prediction for the best fit values (bra 2 = 7.9 x 10~ 5 eV 2 and 
tan 2 (9 = 0.40). 



parameter 




signal variation 


tan 2 9 = 0.26 




+ 13.5% 


tan 2 6 = 0.67 




-12.2% 


bra 2 = 6.9 x 10~ 5 


eV 2 


+0.11 TNU 


bra 2 = 9.3 x 10~ 5 


eV 2 


-0.10 TNU 



as 0.2% at la in Ref. |l2j) translates into a 3a uncertainty of 0.6% on the signal. 



VI. CONCLUDING REMARKS 



We summarize here the main points of this paper: 

(1) Based on a detailed geochemical and geophysical study of the region near Kamioka, 
we have determined the regional contribution to the signal from Uranium-geoneutrinos: 

S TCg = (15.41 ± 2.98) TNU . (20) 

(2) By using global mass balance arguments, we have determined the contribution from 
the rest of the world: 

S RW = [(2.25 + 14.76 X m) ± (-0.55 + 2.61 x m)] TNU , (21) 

where m is the Uranium mass in the Earth, in units of 10 17 kg. 

(3) Our prediction for the signal as a function of m is presented in Fig.Q which shows the 
potential of geo-neutrinos for determining how much Uranium is in the Earth (As discussed 
in the paper, the range of the prediction for a given mass is mostly due to experimental 
determinations of local abundances and to the geometrical distributions of trace elements 
in the mantle). 

(4) Measurements of the antineutrino signal from Uranium can provide crucial tests for 
models. In particular, estimates by different authors for the Uranium mass within the 
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important paradigm of the Bulk Silicate Earth are all in the range (0.7,0.9) x 10 17 kg, which 
translates into a signal 23 < S(XJ) < 31 TNU. 

(5) A full exploitation of the geo-neutrino signal demands that the mixing angle is deter- 
mined more precisely. 
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